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The  research  progratrf  on  "Transport  of  Nonneutral  Plasmas"  includes  both  on¬ 
going  experiments  and  theory  on  transport  in  pure  electron  plasmas;  and  the  design, 
construction  and  operation  of  a  new  pure  ion  plasma  apparatus.  Both  aspects  of  this 
program  are  proceeding  well,  as  detailed  below  and  in  the  scientific  publications 
referenced  herein. 

Recent  experiments  on  the  electron  apparatus  have  elucidated  the  transport  and 
turbulence  associated  with  diocotron  instabilities  and  vortex  structures,  and  have 
measured  the  vortex  pairing  instability  in  detail.  Recent  theory  work  has  analyzed  the 
newly  observed  /  =  1  diocotron  instability,  analyzed  "^negative  temperanire"  2D  thermal 
equilibria,  and  investigated  the  crystalline  equilibria  of  cold  ion  clusters.  The 
characteristics  of  the  new  ion  apparatus  have  been  determined,  the  major  components 
have  been  designed,  and  purchasing  and  construction  is  underway.  We  expect  to  have 
the  new  apparatus  operational  as  planned  near  the  end  of  the  second  year  of  this  grant. 

We  have  been  studying  the  transport  which  results  from  the  essentially  two- 
dimensional  ExB  drift  of  plasma  across  field  lines.  These  motions  are  the  fastest 
radial  transport  process  that  can  take  place  in  our  devices.  The  electron  dynamics  can 
be  approximated  by  2-dimensional  guiding  center  theory.  Here,  the  axial  bouncing  of 
the  individual  electrons  averages  over  any  z -variations  at  a  rate  fast  compared  to  r-0 
motions.  In  this  approximation,  the  2D  drift-Poisson  equations  for  the  evolution  of  the 
electron  column  are  isomorphic  to  the  2D  Euler  equations  for  an  inviscid  fluid  of 
uniform  density  p.^  In  the  electron  system,  the  flow  vorticity  is  proportional  to  the 
density  n ,  which  we  measure  directly. 

Thus,  in  this  approximation  an  initial  distribution  of  electrons  n  (r  ,0,r  =0)  in  a 

cylinder  with  vorticity  Q  cc  n  will  evolve  exactly  the  same  as  an  initial  distribution  of 
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vorticity  Q  in  a  uniform  fluid  such  as  water.  The  electron  experiments  have  many 
advantages  for  testing  2D  fluid  theory:  the  magnetic  field  tends  to  keep  the  electron 
system  2D,  and  the  electron  column  has  low  internal  viscosity  and  has  no  boundary 
layers. 

Diocotron  Instabilities 

When  the  charge  density  profile  nQ(r)  is  non-monotonic  (i.e.  at  least  partially 
"hollow"),  unstable  diocotron  modes  may  appear.  These  modes  vary  as 
5n (r ) exp  ( ikz  +ilQ-icot  } ,  with  ^  =  0  for  the  2D  modes  considered  here. 
Experimentally,  the  unstable  modes  are  found  to  be  distinct  from  the  stable  modes.^’^ 
That  is,  for  any  given  /,  the  stable  and  unstable  modes  coexist,  with  distinct 
frequencies  and  radial  eigenfunctions.  We  observe  instabilities  with  I  >2,  as  expected. 
Furthermore,  we  observe  a  robust  exponential  instability  for  /  =  1  where  none  is 
predicted.  These  instabilities  give  rise  to  large  0-asymmetries  in  n(r,0)  and  0(r,0), 
which  then  result  in  radial  transport  on  tfie  drift  motion  time  scale.  Our  data  for 
nir,Q,t)  characterizes  the  k  =  0  modes  rather  completely. 

For  /  =  1 ,  we  are  able  to  observe  exponential  mode  growth  over  a  range  of  up  to 
2  decades  before  nonlinear  saturation  occurs.^’^’^^^  The  instability  appears  to  be  a 
linear  process,  i.e.  independent  of  the  perturbation  amplitude  at  r=0.  Typically,  the 
instability  e-folds  in  about  6  rotation  periods  for  the  column. 

Prior  theory  had  concluded  that  there  are  no  instabilities  for  /  =  1,  whereas  recent 
theory  work  has  elucidated  an  instability  which  grows  algebraically  with  time.®*^®  The 
theoretical  growth  with  time  can  be  obtained  by  numerically  integrating  a  linear 
Laplace  transform  solution  which  starts  from  initial  conditions  corresponding  to  the 
experiments.  Alternately,  one  can  perform  an  asymptotic  analysis  on  the  Laplace 
transform  solution,  obtaining  a  perturbadon  proportional  to  as  r  -»«>.  Although  the 
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differences  between  experiment  and  theory  here  are  large,  there  are  intriguing 
similarities.  The  early  time  instability  growth  is  seen  to  be  similar.  Further,  the  time- 
asymptotic  theory  perturbation  is  self-shielding,^  as  is  observed  experimentally  at  all 
times.^  These  experimental  and  theoretical  results  on  the  new  /  =  1  instability  on 
hollow  columns  have  recently  been  published  as  two  Physical  Review  Letters, and  as 
contributed  talks  at  APS/DPP-89,^  and  at  APS/DFD-89.^^’^^ 

For  I  =  2,  an  exponentially  unstable  mode  is  observed  for  a  wide  range  of  hollow 
profiles,^’^  and  these  results  are  basically  consistent  with  2D  theory  predictions. 
Numerical  solution  of  the  eigenvalue  equation  for  realistic  smooth  profiles  has  given 
fair  agreement  with  unstable  /  =  2  growth  rates  in  these  preliminary  experiments.  The 
/  =  2  perturbation  grows  exponentially  until  the  mode  is  fully  saturated  and  forms  two 
vortex  structures  in  a  background  of  lower  vorticity.  The  two  vortices  persist  for 
many  rotations,  then  merge  to  the  center.  Presumably,  the  two-vortex  state  is  unstable, 
and  the  merger  is  the  result  of  a  vortex  pairing  instability.  There  may  be  extended 
filamentary  structures  arising  during  the  vortex  formation  and  disintegration,  giving 
rise  to  density  fluctuations  on  smaller  spatial  scales. 

This  work  in  instability-driven  transport  has  been  the  subject  of  several 
invited^"^’^*  and  contributed  talks,^®  and  has  been  published  in  the  conference 
proceedings  of  the  lAEA^  and  a  US-Japan  workshop.^ 

Vortex  Dynamics 

The  two  vortex  merger  process  can  be  studied  as  an  isolated  process,  independent 
of  any  diocotron  instability/  process.  In  a  sense,  this  is  a  generalization  of  our  recent 
work  on  the  (stable)  nonlinear  /  =  1  diocotron  mode,^  which  may  be  thought  of  as  a 
sing!?  vortpx  interacting  with  a  wall  Straightforward  maiiipulaiion  techniques  allow 
us  to  form  an  initial  condition  consisting  of  two  electron  columns  of  chosen  profile  and 
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placement.  Here,  we  consider  the  particular  case  of  two  equal  columns  which  are 
placed  symmetrically  on  either  side  of  the  cylindrical  axis,  separated  by  a  distance  2D. 

We  observe  that  the  behavior  of  the  two  vortices  depends  dramatically  on  their 
separation  to  diameter  ratio.^^’^  If  the  vortices  are  separated  by  more  than  2.0 
diameters,  we  observe  that  they  orbit  around  each  other  relatively  unperturbed  for  up 
to  10“^  orbits.  If  the  vortices  are  initially  separated  by  1.9  diameters,  their  mutual 
interaction  quickly  results  in  filamentary  tail  formation,  but  the  vortices  still  orbit 
around  each  other  for  about  100  orbits  before  merging  at  the  center.  If  they  are 
separated  by  1.8  diameters,  we  observe  merger  at  the  center  in  less  than  one  orbit 
period.  Thus,  the  time  to  merge  abruptly  increases  from  10  |isec  to  about  1  sec  as  the 
separation  varies  from  1.8  to  2.0  diameters. 

This  may  represent  the  cleanest  experimental  measurement  to  date  of  this 
fundamental  vortex  interaction  process.  The  low  inherent  viscosity  of  the  plasma  and 
the  total  absence  of  boundary  layers  at  the  wall  both  contribute  to  absence  of 
secondary  effects  which  could  obscure  its  result.  These  preliminary  results  are  in  fair 
agreement  with  theory  and  computational  results,  and  with  an  experiment  in  a  water 
tank  at  Australian  National  University. 

The  orbit  frequency  and  vortex  shape  distortion  can  also  be  measured  as  a 
function  of  separation  of  the  two  vortices.  We  find  that  the  orbit  frequency  is  well- 
modelled  by  the  interaction  of  2  "point"  vortices  of  appropriate  total  circulation  (i.e. 
total  charge  per  unit  length  axially),  but  that  the  wall  interaction  must  be  included  for 
large  separations.  We  also  find  that  the  interacting  vortices  elongate  towards  each 
other  even  in  the  absence  of  merger.  The  elongation  appears  roughly  consistent  with 
various  "moment  models"  and  numerical  calculations,  but  the  experimental  data  has 
not  yet  been  analyzed  in  detail. 
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Of  course,  2  vortices  may  merge  even  if  they  are  not  symmetric  in  size  or 
placement.  Experimentally,  the  asymmetry  is  easily  controlled.  Also,  it  is  easy  to  take 
3  vortices  as  the  initial  condition.  Here,  the  parameter  space  is  substantially  larger, 
and  chaotic  behavior  is  often  observed. 

This  work  on  vortex  dynamics  has  been  presented  at  Plasma  Physics^  ^  and  Fluids 
Dynamics  conferences,^  and  is  being  prepared  for  publication. 

Negative  Temperature  Equilibria 

Theoretically,  one  candidate  for  the  configuration  resulting  from  vortex  mergers 
or  violent  shear-flow  instabilities  is  the  "most  probable  state"  obtained  by  maximizing 
the  Vlasov  entropy  S  =  -jn  In/i  of  the  distribution  of  guiding  centers,  subject  to  the 
constraints  imposed  by  certain  conservation  laws.  In  the  continuum  limit,  the  integral 
of  any  function  of  density  is  conserved.  Nonideal  dissipation  or  coarse-graining 
eliminates  the  small  scales  of  turbulence,  and  breaks  the  conservation  of  all  but  the 
simplest  of  these  functionals.  Hence  only  a  few  gross  conservation  laws  may  survive. 

We  have  analytically  studied  the  equilibria  which  result  from  taking  the  total 
charge,  the  canonical  angular  momentum,  and  the  electrostatic  energy  to  be  the  only 
conserved  quantities.^’*’^’^  (This  assumes  that  the  guiding-center  degrees  of  freedom 
are  decoupled  from  velocity  space  over  the  time  scales  which  characterize  turbulent 
relaxation.)  The  constrained  maximum-entropy  configurations  are  then  guiding-center 
thermal  equilibria,  described  by  the  equation 

/i(r,0)  = /ioexp{-P(e<l)-a)r2))  ,  (1) 

where  the  parameters  Kq,  P,  and  co  are  determined  by  the  charge,  energy,  and  angular 


momentum. 
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We  find  that  "negative  temperature"  solutions  exist  for  certain  parameter  regimes 
and  that  these  equilibria  are  displaced  from  the  geometrical  center  of  the  system  at 
large  energies  or  angular  momenta.  Such  an  off-axis  state  can  be  regarded  as  a 
dynamic  equilibrium  or  equivalently  as  a  finite-amplitude  diocotron  mode.  1/(3  is  an 
effective  temperature,  and  is  negative  for  the  off-axis  equilibria.  These  results  have 
been  confirmed  by  Monte  Carlo  simulations  of  ensembles  of  guiding  centers.  We  find 
that  the  statistical  fluctuations  can  be  surprisingly  large,  much  as  in  a  system 
undergoing  a  phase  transition.  However  the  relation  between  these  thermal  equilibria 
and  the  final  states  observed  in  experiments  is  not  yet  clear. 

These  off-axis  equilibrium  results  have  been  discussed  at  conferences,^’^*  and 
have  been  published  in  Phys.  Rev.  Lett.^  with  a  more  complete  description  submitted 
to  Phys.  Fluids.^ 

Ion  Crystals 

When  a  single  species  plasma  in  thermal  equilibrium  is  cooled  to  a  sufficiently 
low  temperature,  it  passes  first  into  the  liquid  state  and  then  into  a  crystalline  state. 
The  theory  of  these  states  is  well  understood  for  infinite  sized  systems  with  no 
boundaries:  this  system  is  called  the  One  Component  Plasma  (OCP).  Computer 
simulations  and  analytic  theory  for  an  infinite  homogeneous  OCP  predict  that  for  r>2 
the  system  of  charges  begins  to  exhibit  local  order  characteristic  of  a  liquid,  and  for 
r  >  172  there  is  a  first-order  phase  transition  to  a  bcc  crystal. 

However,  the  physics  is  substantially  more  complicated  for  realistic  clusters  of 
10^-10'*  particles.  In  cryogenic  electron  plasma  experiments,  r-1-2  have  been 
achieved^;  in  ion  plasma  experiments  at  NIST,  T  values  in  the  range  of  several 
hundred  have  been  measured,  putting  the  system  well  into  the  regime  of  strong 
correlation. 
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We  have  studied  the  properties  of  these  finite-sized  clusters  using  a  combination 
of  analytic  theory  and  computer  simulation.'*’^ For  moderate  values  of  F 
(2<r<l(X))  we  find  that  the  density  of  the  cloud  exhibits  spatial  oscillations.  These 
oscillations  have  maximum  amplitude  at  the  cloud  surface  and  decay  back  to  the 
background  density  aiq  increasing  distance  from  the  surface.  The  oscillations  are 

evidence  of  local  order,  and  the  damping  length  is  a  measure  of  the  correlation  length. 
For  r>140,  the  cloud  separates  into  concentric  spheroidal  shells.  In  this  regime 
charges  rarely  move  from  shell  to  shell,  but  they  diffuse  freely  within  the  shell 
surfaces.  Thus,  the  system  might  be  characterized  as  a  liquid  within  the  shells,  and  as 
a  solid  in  the  direction  perpendicular  to  the  shells,  as  in  a  smectic  liquid  crystal.  For 
still  larger  values  of  F  particle  diffusion  within  the  shells  also  approaches  zero  and  an 
imperfect  2-D  hexagonal  crystal  is  formed  within  each  shell. 

We  find  that  the  lattice  structure  in  the  crystalline  regime  is  quite  different  fi’om 
the  bcc  lattice  predicted  for  an  infinite  homogeneous  OCP.  Simulations  have  been 
performed  on  up  to  2000  ions,  and  no  evidence  of  bcc  structure  has  been  observed. 
Recent  theoretical  work  based  on  the  extremely  oblate  "slab"  limit  of  the  bounded 
crystal  predicts  that  the  system  may  have  to  be  quite  large  before  the  minimum  energy 
state  displays  bcc  symmetry,  perhaps  requiring  as  many  as  60  shells  before  the  infinite 
homogeneous  structure  is  recovered.  Thus,  the  addition  of  a  surface  term  to  the  free 
energy  can  change  the  thermal  equilibrium  lattice  structure,  even  if  this  surface  term  is 
relatively  small  compared  to  the  bulk  term. 

Recently,  the  free  energy  of  the  classical  one-component  plasma  has  been 
calculated  analytically  in  the  crystalline  phase  for  both  fee  and  bcc  lattices  to  0(T^), 
where  T  is  the  temperature.'*  By  application  of  thermodynamic  perturbation  theory,  we 
can  explicitly  evaluate  the  effect  of  three  and  four  phonon  interactions  on  the  partition 
function.  Periodic  boundary  conditions  are  applied  to  make  contact  with  previous 
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numerical  work,  in  which  the  O  (T^)  term  was  assumed  to  be  negligible.  We  find  that 
this  term  is  much  larger  than  previously  diought.  This  increases  the  thermodynamic 
stability  of  the  crystal  phase  over  previous  estimates. 

These  results  have  been  the  subject  of  several  invited^^’^^’^^  and  contributed 
talks, and  have  been  published  in  Physical  Review  and  two  Conference 

Proceedings. 

Ion  Apparatus  Development 

We  are  presently  developing  a  new  apparatus  to  contain  a  nonneutral  ion  plasma 
suitable  for  laser  induced  fluorescence  and  optical  tagging  diagnostics.  The  goal  is  to 
obtain  detailed  understanding  of  the  transport  of  particles  across  magnetic  field  lines  in 
various  parameter  regimes.  The  characteristics  of  the  new  apparatus  have  now  been 
determined,  laboratoiy  space  is  being  renovated  for  this  use,  bids  have  been  obtained 
on  the  major  components,  and  other  mechanical  and  electronic  components  are  being 
designed  and  constructed.  At  present,  we  foresee  no  insurmountable  difficulties  in  the 
development  of  this  apparatus.  The  main  uncertainty  at  present  is  in  acquisition  of  the 
laser  system:  the  ONR  did  not  provide  funds  for  this  purpose,  but  we  hope  to  make  a 
start  on  the  laser  diagnostics  using  cost  savings  from  other  categories. 

Laser  induced  fluorescence  (LIF)  has  been  successfully  used  in  plasma  physics 
for  several  years.  The  use  of  this  technique  has  produced  scientific  results  on  the 
fundamental  plasma  physics  of  gas  discharges  and  alkaline  plasmas  and  also  in 
Tokamak  research.  However,  we  are  in  a  parameter  region  significantiy  different  than 
the  existing  experiments.  Our  ions  will  be  in  the  same  temperature  range  as  that  of 
gas  discharge  or  a  Q-plasma,  but  the  density  is  about  two  orders  of  magnitude  smaller, 
leading  to  a  smaller  number  of  detected  photons.  We  intend  to  use  Magnesium  24  as 
the  ion,  with  plasma  parameters  of  >  lO^cm"^  and  T,-  =  0.1  - 10  eV. 
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In  designing  the  LIF  system,  we  have  considered  a  number  of  effects  which 
complicate  the  signal,  including  Doppler  broadening,  power  broadening,  Zeeman 
splitting,  isotopic  shifts,  hyperfine  structure,  and  transit  time  effects.  While  restrictive, 
none  of  these  complications  present  insurmountable  problems. 

For  optical  tagging,  we  need  to  optically  pump  ions  into  a  state  which  has  a  long 
lifetime  compared  to  the  physical  time  scale  of  the  investigated  process.  A  search 
laser  is  then  used  to  detect  the  test  particles  at  different  spatial  locations  in  the  plasma. 
In  addition  to  spatial  information,  information  about  velocity  space  is  also  obtained  if 
the  laser  band- width  is  narrower  than  the  Doppler  band- width.  Two  kinds  of  long 
lifetime  states  have  been  used  previously:  electronic  metastable  states  and  electronic 
spin  states.  Since  the  lifetime  of  the  metastable  ion  state  is  generally  limited  to  several 
tens  of  seconds,  we  are  going  to  use  the  electronic  spin  state.  Since  these  spin  states 
have  a  lifetime  which  is  probably  limited  only  by  spin  exchange,  the  lifetime  should 
be  sufficient  for  our  purposes  when  the  ion  temperature  is  smaller  than  1  eV.  We 
calculate  that  we  will  need  typically  a  90  [is  pulse  with  a  1  mW  CW  single  mode  laser 
to  flip  all  the  spins  along  the  laser  beam  path  inside  of  the  plasma. 

Magnesium-24  was  selected  for  the  following  properties:  1)  an  optical  transition 
from  the  ground  state  in  the  visible  region  or  near  UV;  2)  no  long  lifetime  metastable 
states  which  could  quench  the  fluorescent  signal;  3)  no  hyperfine  structure  which 
would  complicate  the  LIF  signal;  4)  a  non-integer  total  electron  spin  allowing  us  to 
make  use  of  spin  state  tagging  technique.  For  ^'*Mg,  the  interesting  transitions  (e.g. 
3s^S  1^  ->  3p^P  1^)  are  at  about  280  nm. 

With  this  choice,  the  anticipated  plasma  parameters  are  as  follows: 

Magnesium  plasma,  singly  ionized 
Magnetic  field  B  ~  6T 

Density  n  >  10^  cm“^ 

Brillouin  limit  ng  ~  4x  10^  cm~^ 
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lon  temperatures 

Plasma  diameter 
Plasma  length 
Confinement  time 
Debye  length  (1  eV) 
Cyclotron  frequency 
Plasma  frequency 
Rotation  frequency 
Bounce  frequency 
Ion  collision  frequency  (90°) 
Base  pressure 
Ion  neutral  collision  (90°) 


Ti  -O.l^lOeV 
^'^Doppler  5.3  — >  53.  GHz 
IRp  ~  2.  cm 
Lp  -  25.  — >  50.  cm 
X  -  several  hours 
Xq  -  0.23  cm 
Vci  -  3.8  MHz 
Vpi  >  135  kHz 
V,  >  2.4  kHz 
Vjj  <  4  kHz 
ii  -  0.5->10  Hz 
p  <  10"‘°  Ton- 
<  10^  Hz 


We  need  an  intense  tunable  narrow  band  light  source.  It  appears  that  the  only 
possibility  is  a  dye  laser,  since  the  only  other  sufficiently  tunable  laser  is  solid  state 
Ti-Sapphire  at  700  nm  — >  980  nm.  No  existing  dye  works  at  280  nm,  so  we  will  use 
a  dye  at  560  nm  and  frequency  double.  Our  requirements  for  the  laser  are  the 
following: 

Tunable  and  stabilized; 

Computer  controlled  with  capability  for  measuring  the  laser  wavelength; 

Frequency  doubled; 

Band-width  much  smaller  than  the  Doppler  broadening; 

Laser  power  "sufficient." 

We  believe  that  these  requirements  are  best  m.ct  by  a  CW  laser  with  external  frequency 
doubling,  giving  a  power  of  1  mW  at  280  nm. 

The  characteristics  of  the  optical  transport  and  fluorescence  detection  system  have 
been  analyzed.  Considering  the  solid  angle  of  collection,  grid  and  fiber  optic  losses, 
vacuum  window,  filter  and  lens  losses,  and  photomultiplier  tube  quantum  efficiency, 
we  expect  an  overall  efficiency  of  light  collection  around  1.5  x  10“^.  For  our  plasma 
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parameters,  this  results  in  signal-to-noise  ratios  of  10  or  better  for  a  typical  experiment 
when  integration  times  of  1  msec  are  used. 

The  room- temperature  bore  superconducting  magnet  is  perhaps  the  most  basic 
(and  expensive)  single  part  of  the  experiment.  However,  the  characteristics  of  the 
magnet  could  not  be  specified  until  the  laser  diagnostic  characteristics  were  understood 
and  the  test  ion  selected.  For  example,  the  mass  of  the  ion  determines  the  magnetic 
field  required  for  a  given  ion  density.  Similarly,  the  size  of  the  magnet  is  determined 
by  the  ion  Larmor  radius,  the  ion  Debye  length,  the  optical  access  requirements,  etc. 
With  these  questions  answered,  we  have  specified  the  magnet  and  received  bids  from 
vendors.  Our  specified  magnet  has  a  bore  of  8  3/4",  a  uniform  field  length  of  60  cm, 
and  a  field  strength  of  6  Tesla.  We  expect  delivery  of  the  complete  magnet  system  by 
December  1990. 

The  laboratory  layout  was  dependent  on  the  specification  of  the  magnet,  since 
there  are  numerous  safety  and  practical  concerns  with  a  magnet  of  this  strength  and 
size.  We  have  designed  a  workable  laboratory,  and  renovations  funded  by  UCSD 
should  be  completed  by  August  1990. 

At  present,  we  intend  to  obtain  the  ^'*Mg'^  ions  from  a  vacuum  arc  source.  The 
plasma  produced  by  the  source  will  probably  have  a  substantial  level  of  impurities,  an 
unacceptably  high  temperature,  density  jitters  of  about  20%,  and  the  wrong  radial 
profile.  We  intend  to  get  rid  of  the  unwanted  ions  by  heating  with  RF  at  their 
cyclotron  frequency,  increasing  their  larmor  radius  until  they  hit  the  wall.  We  intend 
to  cool  the  ion  plasma  by  collisions  with  neutral  gas  such  as  He.  If  a  fast  magnesium 
ion  has  a  head-on  collision  with  a  helium  neutral  which  is  at  rest,  it  loses  25/49  of  its 
energy.  Thus  in  10  collisions  a  250  eV  plasma  would  cool  to  room  temperature.  This 
requires  about  10  seconds  at  a  pressure  of  2x  10“^  Torr.  We  have  performed  tests  on 
the  effects  of  puffed  gas  loads  on  UHV  pumps,  and  feel  that  this  neutral  load  is 
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acceptable,  especially  if  the  pump  is  valved  off  during  the  gas  puff.  Finally,  we  intend 
to  manipulate  the  plasma  density  profile  and  total  particle  number  by  well-established 
techniques  such  as  "electric  tilt"  perturbations  and  selective  axial  release  of  particles. 
These  various  manipulations  will  be  done  with  computer  sensing  and  control.  This  is 
a  straightforward  extension  of  techniques  in  use  on  present  apparatuses.  Although 
there  may  be  various  difficulties  in  the  development  of  this  new  experiment,  we  are 
optimistic  that  the  new  ion  machine  will  be  in  operation  on  schedule. 
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